Competing Ferromagnetic and Charge-Ordered States in Models for Manganites: 

the Origin of the CMR Effect 



Cengiz §en,^ Gonzalo Alvarez,^ and Elbio Dagotto^'^ 

^National High Magnetic Field Laboratory and Department of Physics, Florida State University, Tallahassee, FL 32310 
'^Computer Science & Mathematics Division and Center for Nanophase Materials Sciences, 
Oak Ridge National Laboratory, Oak Ridge, TN 37831 
'^Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, TN 32831 
^'Department of Physics and Astronomy, The University of Tennessee, Knoxville, TN 37996 

(Dated: February 5, 2008) 

The one-orbital model for manganites with cooperative phonons and superexchange coupling 
Jaf has been investigated via large-scale Monte Carlo (MC) simulations. Results for two-orbitals 
are also briefly discussed. Focusing on electronic density n—0.75, a regime of competition between 
ferromagnetic (FM) metallic and charge-ordered (CO) insulating states was identified. In the vicinity 
of the associated bicritical point, colossal magnetoresistance (CMR) effects were observed. The CMR 
is associated with the development of short-distance correlations among polarons, above the spin 
ordering temperatures, resembling the charge arrangement of the low-temperature CO state. 
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Introduction: The colossal magnetoresistance of the 
manganites is an example of the complex behavior and 
nonlinearities that can emerge in materials where sev- 
eral degrees of freedom are simultaneously active P, Q . 
Understanding theoretically the CMR effect is impor- 
tant not only for its intrinsic value in the Mn oxide con- 
text, but also to provide a paradigm for rationahzing a 
plethora of related complex phenomena unveiled in sev- 
eral other transition metal oxides in recent years 0] . The 
current scenarios to understand the CMR behavior rely 
on the existence of competing states and the emergence 
of nanometer-scale electronic structures 0] . Using simple 
models of phase competition (valid at large length scales) 
and resistor-network approximations, CMR behavior was 
observed Q- However, for a deeper understanding, it is 
imperative to unveil the CMR effect in more fundamen- 
tal microscopic models, using fully unbiased many-body 
techniques to handle the strong interactions. Along this 
line, considerable progress was recently made with the 
report of large magnetoresistance effects in Mn-oxides 
models in the case where the double-exchange (DE) in- 
duced FM metal competes with a FM insulator made 
out of randomly localized polarons @, @i0]- However, a 
complete rationalization of CMR physics needs a CO an- 
tiferromagnetic (AF) state as the true competitor of the 
FM metal, since it is under these circumstances that the 
magnetoresistance effect is truly colossal P, 0] . For this 
goal, it is crucial to incorporate the t2g AF superexchange 
Jaf, known to be fundamental to reproduce the CO-AF 
states observed in real manganite phase diagrams [2]. 

The main purpose of this manuscript is to present 
large-scale computational studies of models for mangan- 
ites, incorporating now explicitly the Jaf interaction. 
The main result is that a CMR effect, with CMR ra- 
tios as large as 10,000%, is unveiled using realistic models 
and electronic densities, and unbiased methods. The Jaf 



coupling is shown to be crucial for the magnitude of the 
effect. A simple picture for the origin of the CMR effect 
is discussed, that relies on nanometer-scale short-range 
order above the Curie temperature Tc, in qualitative 
agreement with previously proposed mixed-phase scenar- 
ios. The comprehensive investigation reported here was 
possible by using hundreds of nodes of the Cray XT3 su- 
percomputer operated by the National Center for Com- 
putational Studies, at Oak Ridge National Laboratory 
(ORNL). The same effort on conventional PC clusters 
would have demanded several years for completion. 

Models and Techniques: The models and methodology 
used in this effort have been extensively discussed be- 
fore Then, only a brief description is here included. 
The one-orbital Hamiltonian is given by H—Huj+Hap, 
where Jfib contains the standard fermionic nearest- 
neighbor (NN) hopping term at infinite Hund's coupling, 
plus the interaction with the oxygen phonons which are 
assumed classical. The full iJib term, includingon-site 
disorder, is explicitly defined in Eq.(2) of Ref. |7|. The 
two-orbital model, also used here, is identical to Eq.(4) 
of Ref. 0- The electron-phonon coupling is denoted by 
A, the strength of the on-site quenched disorder is A, 
the chemical potential /i regulates the density n, and 
the hopping t—1 is the energy unit [1]. In this study, 
the second term Haf=Jaf ' pl^^ys a key role. 
This interaction represents the NN superexchange, with 
coupling Jaf, between the (assumed classical) t2g spins 
Si. The study of this model involves the standard exact 
diagonalization of the fermionic sector for a given classi- 
cal phononic and t2g spins configuration, updated via a 
MC procedure 0]- The charge-charge correlations were 
evaluated using C{j)={l/N) ~ (")^) "^ith 
the electronic density at site i, and N the total number 
of sites. The resistivity p has been calculated by tak- 
ing the inverse of the mean conductivity a, where the 
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latter is related to the conductance G by G=aL'^~'^, d 
being the dimension, and L the linear size of the lattice. 
The conductance G was obtained using the Kubo for- 
mula [i^. For further details, the reader should consult 
Sec. II of Ref. 17|. Finally, note that finite clusters do not 
lead to true singularities, such as a critical temperatures. 
However, decades of investigations in a wide variety of 
contexts have shown that finite systems can accurately 
capture the rapid growth of correlation lengths in nar- 
row temperature ranges expected from criticality. These 
temperatures can be safely associated with true critical 
temperatures, and this is the convention followed here 
and in many previous studies [3, 0, H, 0, 0] . 

Clean-Limit Phase Diagram: The key novelty of this 
investigation is the existence of an insulating CO/AF 
state competing with the FM metal to generate the CMR. 
In Fig. [U the clean-limit A=0 one Cg-orbital model phase 
diagram is presented, at an experimentally realistic den- 
sity n=0.75 and intermediate A=1.2. At Jaf=0, the 
ground state is a DE generated standard FM metal, with 
charge uniformly distributed. With increasing Jaf, a- 
first-order transition occurs to a CO/AF state schemat- 
ically shown in Fig. [1] The holes are arranged reg- 
ularly, separated by distances 2 and \/5. This state is 
degenerate with a state rotated in 90°. The spins arc also 
indicated. This type of states are the analog of the more 
realistic CO / AF states that exist in two-orbital Mn-oxide 
models, such as the CE state 0. 

Clean-limit CMR: The most important result of these 
investigations is shown in Fig. [21 where the p vs. tem- 
perature T curves are shown. Consider Fig. [SJa): here p 
vs. T presents the expected insulating behavior at large 
Jaf , and the canonical bad- metal DE form at small Jaf 
(or reducing A). The remarkable result appears in be- 
tween, mainly in the narrow Jaf interval approximately 
between 0.02 and 0.0325. In this regime, p vs. T presents 
a canonical CMR shape, with insulating behavior at large 
T, transforming into a broad high peak upon cooling (log- 
arithm scale used), followed by metallic behavior at low 
T. The Tq is approximately located at the resistivity 
peak. This is in agreement with the experimental phe- 
nomenology of CMR manganites, and with the theoreti- 
cal scenarios 0, 0] based on a CMR emerging from phase 
competition between a FM metal and a CO-AF insula- 
tor. The result is particularly remarkable considering the 
relatively small clusters used in our studies: the origin of 
the CMR effect unveiled here must lie in phenomena oc- 
curring at the nanometer scale, as recently remarked 11 1 . 
The only "price" to pay is the tuning of couplings: the 
CMR shape is obtained in narrow intervals AA^O.l (inset 
Fig. m^a)) and A Jaf~0.02. However, this delicate tun- 
ing is removed by adding quenched disorder, as shown 
before 0, 0, ^'^^ ^^^^ below in the present effort. 

In Figs. [DIb,c), the effects of magnetic fields H are 
shown. In agreement with experiments, the peak in p 
is strongly suppressed by magnetic fields that are small 
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FIG. 1: (Color online) Clean limit MC phase diagram using 
8x8 and 12x12 lattices, at n=0.75 and A=1.2. The AF/CO 
state is schematically shown, with the radius of the circles 
proportional to the electronic density, and arrows represent- 
ing the i2g spins. Charge is uniform in the competing FM 
state. At each temperature, 1x10^ thermalization and 5x10* 
measurement MC steps were carried out for the 8x8 clus- 
ters (and ~7500 and 5000, respectively, for the 12x12 cluster 
points indicated by red stars) . Random starting spin configu- 
rations were used in all simulations, but the results have also 
been checked for convergence using ordered starting configu- 
rations. Tc is defined as the temperature where the spin-spin 
correlations at the maximum cluster distance nearly vanish. 
Similarly, Tap /co is the temperature where both AF and FM 
correlations vanish at the maximum allowed distances. At 
Tap/co, the charge structure factor n(q) at q = qco also 
vanishes. Inset: energy vs. Jap at very low T~0, with the 
FM (CO/AF) phase in black (red). Green dots indicate a 
G-ordered AF regime. Here, the spins were frozen to their 
expected states and the phonons were MC relaxed. 



in the natural units t~\. The concomitant CMR ratios 
[{p{Q)- p{H)) / p{H)]xlQQ% are as large as 10,000%. The 
comparison between 8x8 and 12x12 clusters also show 
that size effects are mild in these investigations. 

Influence of Quenched Disorder, Dimensionality, and 
Number of Orhitals: Fig. [3ja) shows typical p vs. T 
curves in the presence of on-site-energy quenched disor- 
der. As discussed before [1,0], the disorder enhances the 
tendency toward having p vs. T curves with the canoni- 
cal CMR shape, and the fine tuning problem is avoided. 
This is important to rationalize the universality of the 
CMR effect, present in so many Mn oxides that they 
cannot all be fined tuned to the same couplings. 

To complete our investigations, the effects of dimen- 
sionality and orbital number were also addressed. In 
Fig. [31(b) results are shown using a 4x4x4 cluster, as 
a representative of a 3D lattice, in the clean limit. The 
p vs. T curves are very similar to those found in 2D sys- 
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FIG. 3: (Color online) (a) p vs. T in the presence of quenched 
disorder A. Up to ten different disorder realizations were used 
in calculations with quenched disorder. Only small changes 
between configurations were observed. MC steps and starting 
configurations are as in Fig. [T] (b) p vs. T using a 4x4x4 
lattice, parametrized with A, at Jaf=0.03. (c) Two orbitals 
p vs. T results using an 8x8 lattice for Jaf=0.05. In (b) 
and (c), 4000 thermalization and 4000 measurement MC steps 
were used, n=0.75, and the clean limit A=0 was studied. 



FIG. 2: (Color online) Resistivity p vs T curves for various 
parameters: (a) Fixing A=1.2 and varying Jaf- Arrows in- 
dicate Tc's. Results at A=0.8 and A=1.0, with Jaf=0.0, are 
also shown. Inset: results fixing Jaf=0.03 and varying A. 
(b) Effect of magnetic fields (indicated, in t units) on p us- 
ing Jaf=0.0325, on an 8x8 lattice, (c) Same as (b) but for 
Jaf=0.035, on a 12x12 lattice. In (a) and (b), MC steps and 
the starting configurations are the same as in Fig. [T] In (c) , 
7500 thermalization and 5000 measurement steps were used. 
Typical error bars are indicated in (a). 



terns. In Fig. [Sfc), results using a two-orbital model are 
presented, with similar conclusions. Figure [3]^c) and oth- 
ers not shown, suggest that the results are qualitatively 
the same for the one- and two-orbital models [l^ . 

Understanding the CMR Effect: To understand the 
CMR effect found here, it is important to analyze short- 
distance charge correlations above the ordering tempera- 
tures. In particular, the correlations at distance -\/5, de- 
noted C{V5), are important since they are robust in wide 
ranges of temperatures and couplings, they are ubiqui- 
tous in most MC snapshots, and of course they are also 
characteristic of the CO/AF state in Fig. [TJ A clear 
relation between C{^/5) and p is shown in Fig. lUJa), 
where they are plotted together for a representative Jaf- 
In the regime below a new characteristic temperature 
T*~0.15, both quantities nicely track each other, show- 
ing that the increase of short-distance charge correla- 
tions causes the increase in p upon coolin g. This fig- 
ure is conceptually similar to Fig. 4 of Ref. [iSl, where p 
for Lao.yCao.sMnOa was shown together with the neu- 



tron scattering (1/4, 1/4,0) polaron peak that describes 
short-distance polaron correlations. Such a "correlated 
polaron" signal is the analogous of short-range CO in 
our one-band model, a precursor of the CO state sta- 
bilized at low T close in parameter space. Figure HJa) 
shows that the short-distance charge order is reduced at 
low T but it does not disappear: in this regime the FM 
order of the spins, together with some reduction in charge 
correlations, causes the dramatic resistance drop. 

The scale T* where short-range correlations start to 
build up was already predicted based on phenomeno- 
logical investigations it was also found experimen- 
tally [1J|, but had not been reported before using basic 
models for manganites. The resistivity becomes nearly 
flat above T*. In our investigations, a pseudogap (PC) 
in the density-of-states A^(w) was also observed [l^, and 
its inverse at pc is shown in Fig. [DJa). The PC disap- 
pears at a higher temperature T^', probably related to 
polaron formation. The resistivity is not affected by , 
but upon further cooling to T* those polarons become 
correlated and p dramatically increases [l6| . 

Finally, while averaged quantities are certainly crucial 
for quantitative studies, the "by eye" examination of MC 
snapshots is also illustrative. In Fig.Ufb), a typical snap- 
shot in the CMR regime is shown. Pairs of holes located 
at the distances and 2 are highlighted. Clearly, MC 
equilibrated configurations do not have randomly dis- 
tributed holes (i.e., it is not a gas of heavy polarons), 
but special distances are preferred over others: the po- 
larons are correlated. Those distances are the same that 
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FIG. 4: (Color online) (a) MC aver- 
aged C(\/5) vs. r, showing a qual- 
itative similarity with the rescaled 
resistivity (also shown). This agree- 
ment occurs below the T* indi- 
cated. At higher T, p is flat and 
C(\/5) nearly vanishes. Also shown 
is the inverse of N{uj = fj,), to in- 
dicate the PG formation at T^K 
(b) Typical MC snapshot with the 
radius of the circles proportional 
to the local charge density. Also 
shown are the hole-hole distances 
\/5 and 2 of relevance (see text). 



characterize the low-T CO state (Fig. [T]). Also, in the 
MC snapshots the t2g spins at the hole locations and 
their four neighbors were found to be polarized similarly 
to that shown in Fig. [T] Thus, puddles of the CO-AF 
state appear in the MC snapshots, and their existence is 
correlated with the shape of the p vs. T curves. 

Conclusions: The results reported here show that re- 
alistic models for manganites - with the DE, electron- 
phonon, and Jaf interactions - can explain the CMR 
effect and its magnitude. The coupling Jaf is crucial, 
since it is needed to stabilize the CO-AF state that com- 
petes with the FM metal. The origin of the CMR is 
the formation of nano-scale regions above Tq, with the 
same charge and spin pattern as the low-T insulating 
CO-AF state. To obtain CMR effects, clusters of just 
a few lattice spacings in size appear sufficient 
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results show that the CMR effect is much larger when 
the insulating competing state has charge ordering ten- 
dencies, as opposed to having merely randomly localized 
polarons. The latter state has been recently identified 
using STM techniques, in a bilayered a;=0.3 manganite 
with mild in-plane magnetoresistance [l3|. If the same 
STM experiment is repeated at x=OA, where the CMR is 
much stronger, it is predicted that aggregates of polarons 
resembling a CO state should be observed above Tq. 
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